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Abstract. - Non-magnetic impurities break a quantum spin chain into finite segments and induce 
Friedel-like oscillations in the local susceptibility near the edges. The signature of these oscil- 
lations has been observed in Knight shift experiments on the high-temperature superconductor 
YBa2Cu306.5 and on the spin-chain compound Sr2Cu03. Here we analytically calculate NMR 
spectra, compare with the available experimental data for Sr2Cu03, and show that the inter- 
chain coupling is responsible for the complicated and so far unexplained lineshape. Our results 
are based on a parameter-free formula for the local susceptibility of a finite spin chain obtained 
by bosonization which is checked by comparing with quantum Monte Carlo and density-matrix 
renormalization group calculations. 



Introduction. — An important tool to study the lo- 
' cal spin dynamics in strongly correlated electron systems 
. is nuclear magnetic resonance (NMR) . NMR experiments 
' have been instrumental in investigating spin fluctuations 
, and impurity effects in high-temperature superconductors 
' [1], as well as in confirming the triplet nature of supercon- 
. ductivity in Sr2Ru04 [2]. Quite recently, NMR was also 
' used to study the CuO chains in YBa2Cu306.5 (YBCO) 
, [3]. The NMR study showed that the chain ends induce 
' Friedel-like oscillations which manifest themselves also in 
, the Cu02 planes. Similar oscillations have also been ob- 
' served earlier in the prototypical quasi one-dimensional 
'S = 1/2 spin chain compound Sr2Cu03 (SCO) [4,5]. 
Theoretically, a large alternating component of the lo- 
cal susceptibility near the end of a semi-infinite Heisen- 
berg chain has been predicted [6]. Other studies (for a 
recent review see ref. [7]) have addressed local spin corre- 
lations near a chain end by numerical means in a variety of 
one- dimensional models ranging from the frustrated and 
dimerized spin-1/2 chain to spin ladders and the spin-1 
Heisenberg chain [8]. 

An interesting open question concerning the physics of 
the spin-1/2 Heisenberg chain is whether its transport 
properties are ballistic or diffusive. The theoretical re- 
sults are contradictory [9-14] but seem to point to ballistic 



transport perhaps related to the integrability of the model 
by Bethe ansatz. NMR and muon spin relaxation exper- 
iments on SCO [15, 16], on the other hand, have found 
diffusive behavior. In order to analyze these experiments 
on a quantitative level in the future, it is first of all im- 
portant to know in how far a spin-only model for SCO is 
valid. In particular, it has been claimed in ref. [5] that 
already the NMR spectra - where only static correlations 
are tested - can only be understood if a couphng to the 
lattice is taken into account. 

In the first part of this letter we will calculate NMR 
spectra for a Heisenberg chain with a Poisson distribu- 
tion of non-magnetic impurities. This is known to be the 
relevant model for SCO, with chain breaks caused by the 
presence of excess oxygen [17, 18]. We will start with the 
ideal chain but will then show that the interchain cou- 
plings are essential to fully explain the experimental data 
for SCO [4,5]. Our main findings are that the theoretically 
calculated spectra for a spin-only model of weakly cou- 
pled Heisenberg chains are in perfect agreement with ex- 
periment whereas the phenomenological model proposed 
in [5], involving some coupling to lattice degrees of free- 
dom, cannot - if taken seriously - explain the data. The 
calculated NMR spectra are based on a parameter-free for- 
mula for the local susceptibility of finite S — 1/2 XX Z 
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chains at finite temperatures obtained by a bosonization 
approach as explained in the second part of this letter. 
The analytical results allow for a full impurity averaging 
which would be impossible to achieve at low temperatures 
by numerical calculations. 
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Fig. 1: (Color online) Theoretical NMR spectra (3) for a 
Poisson distribution with p = 0.3% non-magnetic impuri- 
ties, parameters J, yl", A^^ as appropriate for SCO (see text), 
r = 4 • 10"*, and /iL = 7.61 T (dashed line). The dotted 
lines mark the maxima corresponding to the bulk susceptibil- 
ity. Insets (a) and (b) show local susceptibilities at various 
temperatures indicated on the plot. 



NMR spectra. — The Hamiltonian of the spin-1/2 
XX Z model with N sites and open boundary conditions 
(OBCs) is given by 

N-l N 

j=i i=i 

(1) 

Here J is the exchange constant, A e [0, 1] an exchange 
anisotropy, and h the applied magnetic field. Due to the 
OBCs, translational invariance is broken leading to a po- 
sition dependent local susceptibility 

d 1 

Xj = -Qj^{SJ)h=a = — (5'|S'tot)'i=o (2) 

where T is the temperature and S^^^ = Sj. The 
hyperfine interaction couples nuclear and electron spins. 
For a chain segment of length N this leads to the 
Knight shift of the nuclear resonance frequency K^^'^ = 

{le/ln) X^j^\ where 7e (7n) is the electron (nu- 

clear) gyromagnetic ratio, respectively. The hyperfine in- 
teraction is short ranged so that usually only and A^^ 
matter. The NMR spectrum is proportional to the distri- 
bution of Knight shifts. Let us assume in the following a 
Poisson distribution of non-magnetic impurities with con- 
centration p and a Lorentzian lineshape with width T for 
each Knight shift. The normalized probability distribution 



is then given by 



N=l 



N 



^ 1 
V - 



(3) 



As we will show in the second part of this letter, bosoniza- 
tion allows us to derive a parameter-free result for xf^"^ iii 
the limit T/ J <C 1 and 1. Because the deviations for 
very small chain lengths are not important for the NMR 
spectra as long as the probability of having such tiny seg- 
ments is low, the only parameters entering in (3) are the 
material-dependent constants J, and A^^ . 

NMR measurements have been performed on the 
Heisenberg (A = 1) chain compound SCO [4,5]. Chain 
breaks in this system are believed to be caused by ran- 
domly distributed excess oxygen leading to the formation 
of Zhang- Rice singlets [17, 18]. From measurements of 
the total susceptibility it follows that J 2200 K [17]. 
By a comparison with YBa2Cu306+5 [19] and theory 
the hyperfine couphng constants A°/(2?i7„) w —13 T, 
A0f,/(2ft7„) « 2 T, and AV(2;i7„) « 4 T are obtained. 
Here the index denotes the magnetic field direction. We 
calculate the spectra as a function of /i = (1 + K)h^^^ 
where h'^^^ — v/^n is the resonance field for an isolated 
^^Cu atom. In experiment ;/ = 86 MHz [4] and 7„ w 11.3 
MHz/T [20] leading to w 7.61 T. Exemplarily, we 
show the evolution of the lineshape for an ideal chain 
with impurity concentration p = 0.3% in fig. 1. At high 
temperatures a central peak dominates whose position 
corresponds to the bulk susceptibility value. In addi- 
tion, broad edges are visible whose separation increases 
~ h^^^^JvjT\Ti'^'^{v /T) with decreasing temperature with 
V being the spin velocity. These edges are caused by the 
extrema in the local susceptibilities of chain segments with 
lengths N TTv/T (see Figs. 1 (a) and (b), respectively). 
Furthermore, we observe a gradual transfer of weight from 
a peak at high temperatures to a peak corresponding to 
zero Knight shift at low temperatures stemming from the 
increasing number of even chain segments with A^ <C ttv/T 
which become frozen into their singlet ground state (see 
fig. 1(b)). The odd chain segments with A^ <C ttv/T, on 
the other hand, will yield large Knight shifts (see fig. 1(a)). 
This leads to a background which grows in intensity and 
expands with decreasing temperature. The temperature 
T* where the peak corresponding to zero Knight and the 
peak corresponding to the bulk susceptibility value have 
equal height is therefore directly related to the impurity 
concentration. For p ^ 1 we analytically find 



^ „ 3-2T*/(,77r= 



(4) 



which is a simple criterion to determine the impurity con- 
centration from the NMR spectra alone. 

The observation of a central peak at high temperatures 
and broad edges with separation Ah ~ h'^^^y/v/T is in 
agreement with experimental observations [4, 5] as shown 
in fig. 2. It has indeed already been pointed out in [4] - 
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Fig. 2: (Color online) (a) Comparison between experimental 
data taken from [4] and theory with material-dependent pa- 
rameters as given in the text and p = 5 x 10"'', F = 4 x 10~* 
at T = 60 K. The peak height has been adjusted to agree with 
experiment and /ij^g = 7.5955. (b) Same as (a) at T = 30 K. 
(c) Calculated oscillations Xj at T = 30 K in an infinite chain 
with an impurity at a; = (chain I) and reflections in the infi- 
nite chain II without impurities due to an interchain coupling 
with strength J± = 5 K. 



based on the theoretical results for the semi-infinite chain 
by Eggert and Affleck [6] - that the edges are a conse- 
quence of the maxima in the local susceptibility. However, 
at temperatures T < 30 K additional structures are visible 
in the experimental spectra. The peak develops shoulders 
whose separation is denoted by Dh in fig. 2(b) following 
the notation introduced in [4]. Furthermore, a splitting of 
the peak, Sh, very different from the weight transfer with 
temperature shown in fig. 1, is observed. In ref. [5] it has 
been tried to explain these features by a phenomenological 
model of mobile bond defects. In particular, the feature 
Sh was ascribed to a periodic arrangement of bond defects 
leading to chain segments of odd length A'' only. However, 
analytically we find that the splitting of the central peak 
would then grow like ~ h°^^J\n/'^{N)/{TN^/^). Due to 
the prefactor ^ J this predicts a very rapid increase of 
the splitting not observed in experiment showing that this 
model is incorrect. Instead, as we will show below, the ad- 
ditional features Dh and Sh are a consequence of the inter- 
chain coupling. The interchain coupling along one of the 
crystal axes perpendicular to the chains is of order J± ^ 5 
K while it is three orders of magnitude smaller along the 
other direction [17,21]. At T ~ 30 K we might therefore al- 
ready expect significant effects of the interchain couplings 
which, however, can be included perturbatively. Doing so 
we find that the susceptibility oscillations near an impu- 
rity residing in one chain (zeroth order) lead to substantial 
reflections in the neighboring chain (first order) for pa- 
rameters appropriate for SCO as shown in fig. 2(c). The 
ratio of the maximum in chain I to the maximum in chain 
II is ~ Ah/{Dh), i.e., the shoulders of the peak, Dh, are 
caused by the maxima of the reflected oscillations. We find 



Dh - h°^^J±y/vT-^/^ but logarithmic corrections can dis- 
guise this scaling as will become clear later on. The split- 
ting of the central peak Sh, see fig. 2(b), has a more com- 
plicated origin. First, there are also a reflections in next- 
nearest neighboring chains (second order). The maxima 
would then yield a splitting Sh ~ /i°gg J^-^/vT"^/^ (again 
ignoring logarithmic corrections). However, for impurity 
concentrations p ^ 5 x 10"** relevant for the experiments 
there is another effect which actually dominates: Includ- 
ing the first and second order reflections from neighbor- 
ing chains, the chain of average length N ^ \/p will not 
have any sites left which show bulk behavior. This means 
that at temperatures T < 30 K the probability for hav- 
ing Knight shifts corresponding to values close to the bulk 
susceptibility starts to decrease dramatically thus leading 
to a drop in intensity in P{K). Therefore Sh denotes not 
a splitting of the central peak but rather a dip of intensity 
at the bulk susceptibility value. The oscillations, which 
now basically spread over the entire crystal, might get fur- 
ther stabilized by anisotropic exchange terms. This might 
also explain the small differences in the lineshape near the 
peak for h \\ a and h\\ c [4,5]. A detailed analysis of these 
anisotropy effects is beyond the scope of this letter. 

We also want to stress that having additional structures 
in the spectra due to interchain couplings is very differ- 
ent from a scenario where instead such structures occur 
due to a direct hyperfine interaction or a dipolar cou- 
pling of the nuclear spins with electron spins in adjacent 
chains. The latter case would lead to Dh/Ah = const. 
{Sh/Ah = const.) as a function of temperature if the 
shoulders (splitting) are caused by this mechanism, re- 
spectively. Only three (two) different temperatures are 
presented in [5] where the shoulders (splitting) are visible, 
respectively, making a detailed analysis impossible. How- 
ever, we find that the data are consistent with Dh/Ah ~ 
1 /T as expected if interchain coupling dominates but cer- 
tainly not consistent with either Dh/Ah = const, or 
Sh/Ah = const.. 

To calculate the NMR lineshapes for SCO at low tem- 
peratures, we have to deal with a two-dimensional array 
of weakly coupled chains. In general, this is a very com- 
plicated task requiring a two-dimensional impurity aver- 
aging. However, for small impurity concentrations signifi- 
cant simplifications are possible. At a given temperature, 
the oscillations extend over a characteristic length ^. If 
N — 1/p ^ ^ then the probability of having two impuri- 
ties in neighboring chains so close to each other that the 
zeroth order oscillations in the chain and the refiections 
from neighboring chains overlap is small. We therefore as- 
sume that reflections in a chain of length N only occur in 
regions where the chain shows bulk behavior. In a chain 
of length N, 2pN reflections from the nearest neighbor 
chains and 2pN reflections from nearest-neighbor chains 
will occur on average. If a chain segment is long enough, 
we consider the zeroth, first, and second order oscillations 
as independent entities. If the segment is too short, we 
reduce the extend of the first and second order oscilla- 
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tions mimicking the overlap. For the chain segments we 
then do the full impurity averaging (3). The theoretically 
calculated lineshapes we obtain this way are in excellent 
agreement with experiment as shown in fig. 2. For both 
temperatures we use the peak to adjust the intensity. For 
T = 60 K we also used p and /i|?gg as parameters and find 
p = 5 X 10-* and h°^^ = 7.5955 T. In fig. 2 (b) we use the 
same values. An even better agreement would be obtained 
here if we choose /i^^g = 7.598 T which is equivalent to a 
deviation of 0.03% from our prediction for the evolution 
of the bulk susceptibility value. Furthermore, the inten- 
sity of the edges Ah is overestimated. This is most likely a 
consequence of our assumption that the zeroth order oscil- 
lations do not overlap with the reflections. Configurations 
where such an overlap occur would wash out the maxima 
of the zeroth order oscillations. In addition, this might 
also point to some deviations from a Poisson distribution 
with short chains occurring less frequently than expected. 



Local susceptibility. 



We now explain how the 



parameter-free results for xf^"^ have been obtained. In 
the low-energy limit, the spin operators can be expressed 
in terms of a boson $ as 



-IV cos v/27rXi,$ 



(5) 



Here is the Luttinger parameter and c the amplitude 
of the alternating part. The integrability of model (1) by 
Bethe ansatz allows it to determine and c exactly for 
all A. Ignoring bulk and boundary irrelevant operators, 
the Hamiltonian (1) is equivalent to a free boson model 



H = 




L+a 

dxd^^{<o) 



where v is the spin velocity, L — Na, and a the lattice 
constant. The bosonic fields obey the standard commu- 
tation rule [<I>(a;),n(a;')] = i6{x - x') with 11 = v'^dt^. 
To calculate the local susceptibility (2) we use a mode 
expansion 



^{x^ja,t) = 




(7) 



sin(.nj7(iV-M))/ 



E 

n=l 



which incorporates the OBCs. Here hn is a bosonic anni- 
hilation operator. Eq. (7) is a discrete version of the mode 
expansions used in [18, 22] with x = ja becoming a con- 
tinuous coordinate for a — s- 0, A'^ — > oo with L = Na fixed. 
Using this mode expansion, the local observables respect 
the discrete lattice symmetry j N + 1 — j corresponding 
to a refiection at the central bond (site) for even (odd), 
respectively. The sites and A^ -|- 1 are added to model 
(1) and we demand that the spin density vanishes at these 
sites. Therefore the upper boundary for the integrals in 
(6) is L -I- a. The zero mode part (first line of eq. (7)) 
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Fig. 3; (Color online) Xj for A = 0.3 and T/J = 0.02. Compar- 
ison between QMC and field theory for iV = 100 and N = 99. 



= iStot ^nd the oscillator 
0, A 1 as 



fulfills ^^.si^y'^/o^ 

part (second line of eq. (7)) vanishes for j 
required. 

Using (5) in the formula for the local susceptibility (2), 
we find Xj = Xj"' + (~l)"'Xj*- The uniform part is inde- 
pendent of position, x""' = X""' = {iStot?}/iTN) and 
we can therefore directly use the parameter-free result de- 
rived in [18]. For the staggered part, on the other hand, we 
find xf = f (cos y2iirr^')osc(cos V27ri?r$S't'ot)zm where 
we have split the correlation function into an oscillator 
and a zero mode part according to (7). Using the cumu- 
lant theorem for bosonic modes (exp (± i-v/27r AT^ $ ) ) osc = 
exp{—TrKL {(/>(/)) osc) we obtain, following [18,23], 



Kl/2 ^^3Kr./2 ^^-fl^ 



(cos ^2t:Kl^)osc = 



(8) 

Here ri{x) is the Dedekind eta- function and Oi(u,q) the 
elliptic theta- function of the first kind. For the zero mode 
part we find 



(cos v^STrATi^^toJ 

zm 

ra sin[27rm_7/(A^ ■ 



(9) 



l)]e 



--nvm^ /{KlLT) 



p-T^vm,-^ /(KlLT) 



with m running over all integers (half-integers) for A^ even 
(odd), respectively. In the thermodynamic limit, A^ — > oo, 
we can simplify our result and obtain 



[;^-nh(-^)] 



\Kj2 



(10) 



with X = ja. This agrees for the isotropic Heisenberg case, 
Kl = 1, with the result in [6]. The amplitude c can be 
determined with the help of the Bethe ansatz along the 
lines of ref. [24]. This leads to c = Az/2 with Az as 
given in eq. (4.3) of [24]. Our result for the staggered 
part of the local susceptibility is therefore parameter free. 
This means that we can directly compare our analytical 
result for Xj with quantum Monte Carlo (QMC) data. For 
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Fig. 4: (Color online) Xi for A = 1 and iV = 100. Com- 
parison between QMC (closed symbols) and field theory (open 
symbols) for T/J = 0.02, 0.05, 0.1, as indicated on the plot. 



an anisotropy A — 0.3, shown in fig. 3, the agreement is 
excellent. 

Next, we come to the experimentally most relevant 
isotropic case, A = 1 {Kl — 1). Umklapp scattering is 
then marginally irrelevant and the scaling dimensions of 
correlation functions have to be replaced by renormaliza- 
tion group improved versions. For Xj the calculations are 
rather similar to those in ref. [25] for the longitudinal spin- 
spin correlation function. We find that we have to replace 
1/Kl 1 — 5 in (9) whereas Kl = 1 in the oscillator part, 
eq. (8). The renormalization of K]^ for this part is incor- 
porated into an effective amplitude c {2TT^g)~^^^ . The 
running coupling constant g depends, in general, on the 
three length scales a;, L, and v/T. At low enough energies 
the smallest scale will always dominate and g is given by 
the solution of 1/g -I- ln{g)/2 = In {Co min[x, L — x,v/T]} 
where Co is a constant. In fig. 4 a comparison between this 
analytic result and QMC data is shown with Co = 6. We 
note that fitting the constant Co improves the results near 
the boundaries. For low temperatures and x, L — x 3> 1, 
however, the value of Cq becomes irrelevant and our re- 
sult for Xj therefore again parameter-free. The agreement 
with QMC is not as good as for the anisotropic case. This 
is a consequence of the fact that g has been derived in the 
limit x,L,v/T ^ 1. However, the deviations are only of 
the order of a few percent and have very little effect on 
the NMR spectra presented in the first part of this letter. 

Finally, the first order refiection xf^^^^ of susceptibility 
oscillations in a neighboring chain, shown in fig. 2(c), is 
given in first order perturbation theory in J± by 



st(l) 



T 



(11) 



Here eq. (10) has to be used for x^ ^^nd G 



1 K 



{S^S^r = cV[;^sinh(2f |j-fc|)]" is the staggered 
part of the bulk two-point correlation function. While 
it would be extremely difficult to obtain accurate numer- 
ical data for two weakly coupled chains at temperatures 
T/J ~ 0.01 as considered in fig. 2(c) we can easily check 
formula (11) at higher temperatures but still T/J ^ 1. 



Fig. 5: (Color online) xf for A = 1, Jx = 0.03J, and T = O.IJ. 
Comparison between DMRG data (closed symbols) and field 
theory (open symbols) for an infinite chain with a non-magnetic 
impurity ai j = (squares) and an infinite neighboring chain 
without impurities (circles). 



Particularly suited to study the case of an infinite chain 
with a single non-magnetic impurity at the origin which 
is weakly coupled to an infinite chain without impurities 
is the density-matrix renormalization group (DMRG) ap- 
plied to transfer matrices. This algorithm allows it to di- 
rectly obtain results in the thermodynamic limit [26,27]. 
In fig. 5, DMRG data are compared to the field theoretical 
formulas (10,11) and good agreement is found. The max- 
imum of Ixj'^^'' I and therefore the separation of the shoul- 
ders Dh scales like T^^/"^ for — 1 with complicated 
logarithmic corrections coming in through the amplitude 
c. This might make it hard to detect this power law in ex- 
periment. The second order reflections can be calculated 
analogously. 

Conclusions. — To conclude, we have derived an ana- 
lytic formula for the local susceptibility of a finite Heisen- 
berg chain. This allows us to calculate NMR spectra for 
spin chains with arbitrary impurity concentrations and 
distributions which would be impossible by numerical cal- 
culations at temperatures T/J ^ 1. We also showed how 
to calculate NMR spectra for weakly coupled spin chains if 
the impurities are dilute. For SCO we have demonstrated 
excellent agreement between our theory and experiment 
showing that SCO is indeed a prototypical quasi one- 
dimensional spin chain compound. More generally speak- 
ing, we have shown that NMR spectra are extremely useful 
to extract information about the impurity concentration 
as well as about the magnetic couplings. In particular, 
the coupling between the chains leads to an additional 
structure and its position allows it to directly extract the 
coupling strength. We also want to remark that this struc- 
ture is very sensitive to the type of interchain coupling. If 
two chains are coupled by a zigzag-interchain coupling, as 
is the case, for example, for SrCu02 [17], Xj*^^^ would be 
zero, i.e., there would not be any reflections to first order 
in neighboring chains. An analysis of NMR spectra can 
therefore also help to clarify the geometry of the relevant 
magnetic exchange couplings. Furthermore, we expect the 
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results presented here to be also helpful for a more detailed 

analysis of NMR spectra for systems like YBCO [3] where 
CuO chains and Cu02 planes are weakly coupled. 

* * * 

The authors thank I. Affleck and, in particular, S. Eg- 
gcrt for helpful discussions about the role of interchain 
couplings. 
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